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The structure of a river network may be seen as a discrete set of nested sub-networks built out 
of individual stream segments. These network components are assigned an integral stream order 
via a hierarchical and discrete ordering method. Exponential relationships, known as Horton's 
laws, between stream order and ensemble-averaged quantities pertaining to network components 
are observed. We extend these observations to incorporate fluctuations and all higher moments 
by developing functional relationships between distributions. The relationships determined are 
drawn from a combination of theoretical analysis, analysis of real river networks including the 
Mississippi, Amazon and Nile, and numerical simulations on a model of directed, random networks. 
Underlying distributions of stream segment lengths are identified as exponential. Combinations of 
these distributions form single-humped distributions with exponential tails, the sums of which are 
in turn shown to give power law distributions of stream lengths. Distributions of basin area and 
stream segment frequency are also addressed. The calculations identify a single length-scale as a 
measure of size fluctuations in network components. This article is the second in a series of three 
addressing the geometry of river networks. 

PACS numbers: 64.60.Ht, 92.40.Fb, 92.40.Gc, 68.70,+w 



I. INTRODUCTION 

Branching networks are an important category of all 
networks with river networks being a paradigmatic exam- 
ple. Probably as much as any other natural phenomena, 
river networks are a rich source of scaling laws [Q, [§ • 
Central quantities such as drainage basin area and stream 
lengths are reported to closely obey power-law statis- 
tics |, §, §, |, §, §, @, §. The origin of this scaling 
has been attributed to a variety of mechanisms includ- 
ing, among others: principles of optimality 0, JJ, self- 
organized criticality BlQ l , in vasion percolation pd| , and 
random fluctuations j^jT 1 2 : |1| One of the difficul- 
ties in establishing any theory is that the reported values 
of scaling exponents show some variation H 7], . 

With this variation in mind, we have in (l£ ] extensive- 
ly examined Hack's law, the scaling relationship between 
basin shape and stream length. Such scaling laws are 
inherently broad-brushed in their descriptive content. In 
an effort to further improve comparisons between theo- 
ry and data and, more importantly, between networks 
themselves, we consider here a generalization of Horton's 
laws JTrJ, [18) . Defined fully in the following section, Hor- 
ton's laws describe how average values of network param- 
eters change with a certain discrete renormalization of 
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the network. The introduction of these laws by Horton 
may be seen as one of many examples that presaged the 
theory of fractal geometry [|l9). In essence, they express 
the relative frequency and size of network components 
such as stream segments and drainage basins. 

Here, we extend Horton's laws to functional relation- 
ships between probability distributions rather than sim- 
ply average values. The recent work of Peckham and 
Gupta was the first to address this natural generaliza- 
tion of Horton's laws [Q. Our work agrees with their 
findings but goes further to characterize the distributions 
and develop theoretical links between the distributions of 
several different parameters. We also present empirical 
studies that reveal underlying scaling functions with a 
focus on fluctuations and further consider deviations due 
to finite-size effects. 

We examine continent-scale networks: the Mississippi, 
Amazon, Congo, Nile and Kansas river basins. As in [16), 
we also examine Scheidegger's model of directed, random 
networks Both real and model networks provide 

important tests and motivations for our generalizations 
of Horton's laws. 

We begin with definitions of stream ordering and Hor- 
ton's laws. Thereafter, the paper is divided into two main 
sections. In Section III, we first sketch the theoretical 



generalization of Horton's laws. Estimates of the Horton 
ratios are carried out in Section IV and these provide 
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basic parameters of the generalized laws. Empirical evi- 
dence from real continent-scale networks is then provid- 
ed along with data from Scheidegger's random network 
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model in Section |v|. In Section VI we derive the high- 
er order moments for stream length distributions and in 
Section VII, we consider deviations from Horton's laws 
for large basins. In the Appendix |a|, we expand on some 
of the connections outlined in Section |y|, presenting a 
number of mathematical considerations on these gener- 
alized Horton distributions. 

This paper is the second in a series of three on the 
geometry of river networks. In the first we address 
issues of scaling and universality and provide further 
motivation for our general investigation. In the third 
article of the series [^lj we extend the work of the present 
paper by examining how the detailed architecture of river 
networks, i.e., how network components fit together. 



II. STREAM ORDERING AND HORTON'S 
LAWS 

Stream ordering was first introduced by Horton in an 
effort to quantify the features of river networks [|f7|. 
The method was later improved by Strahler to give 
the present technique of Horton-Strahler stream order- 
ing p2| . Stream ordering is a method applicable to any 
field where branching, hierarchical networks are impor- 
tant. Indeed, much use of stream ordering has been made 
outside of the context of river networks, a good example 
being the study of venous and arterial blood networks 
in biology ||,|l[^f§[f^|§§§|3. We describe 
two conceptions of the method and then discuss empirical 
laws defined with in the context of stream ordering. 

A network's constituent stream segments are ordered 
by an iterative pruning. An example of stream ordering 
for the Mississippi basin is shown in Figure [|. A source 
stream is defined as a section of stream that runs from 
a channel head to a junction with another stream (for 
an arboreal analogy, think of the leaves of a tree). These 
source streams are classified as the first order stream seg- 
ments of the network. Next, remove these source streams 
and identify the new source streams of the remaining net- 
work. These are the second order stream segments. The 
process is repeated until one stream segment is left of 
order fl. The order of the network is then defined to be 

n. 

Once stream ordering on a network has been done, a 
number of natural quantities arise. These include n u , the 
number of basins (or equivalently stream segments) for 
a given order u>; the average main stream length; 



the average stream segment length; 



the 



average basin area; and the variation in these numbers 
from order to order. Horton |Q and later Schumm 
observed that the following ratios are generally indepen- 
dent of order lu: 
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FIG. 1: Stream segments for u> — 8 up to uj — £1 = 11 for 
the Mississippi River. The spherical coordinates of latitude 
and longitude are used and the scale corresponds to roughly 
2000 km along each axis. 



tions of stream segment lengths l u — Ylv=i we have 
that the Horton ratio for stream segment lengths R,( s ) is 
equivalent to i?/. Because our theory will start with the 
distributions of Zw, we will generally use the ratio R^s) 
in place of Ri . 

Horton's laws have remained something of a mystery 
in geomorphology — the study of earth surface processes 
and form — due to their apparent robustness and hence 
perceived lack of physical (or geological) content. How- 
ever, statements that Horton's laws are "statistically 
inevitable" |3^], while possibly true, have not yet been 
based on reasonable assumptions ||. Furthermore, many 
other scaling laws can be shown to follow in part from 
Horton's laws Thus, Horton's laws being without 
content would imply the same is true for those scaling 
laws that follow from them. Other sufficient assump- 
tions include uniform drainage density (i.e., networks are 
space-filling) and self-affinity of single channels. The lat- 
ter can be expressed as the relation [jjj ^2|, [|| [34| 



I oc Lt 



(2) 



where L» is the longitudinal diameter of a basin. Scaling 
relations may be derived and the set of relevant scaling 
exponents can be reduced to just two: d as given above 
and the ratio In R,( s ) I hii? n |§- Note that one obtains 
R a = R n so that only the two Horton ratios R n and 



(1) Ri(s) are independent. Horton ratios are thus of central 
Since the main stream length averages l u are combina- importance in the full theory of scaling for river networks. 
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III. POSTULATED FORM OF HORTON 
DISTRIBUTIONS 

Horton's laws relate quantities which are indexed by 
a discrete set of numbers, namely the stream orders. 
They also algebraically relate mean quantities such as 
a u . Hence we may consider a generalization to func- 
tional relationships between probability distributions. In 
other words, for stream lengths and drainage areas we 
can explore the relationships between probability distri- 
butions defined for each order. 

Furthermore, as we have noted, Horton's laws can be 
used to derive power laws of continuous variables such as 
the probability distributions of drainage area a and main 
stream length I |7[ [s[ g5): 



P(a) cx a 



and 



p(i) cx r 7 . 



(3) 



These derivations necessarily only give discrete points of 
power laws. In other words, the derivations give points as 
functions of the discrete stream order lo and are uniform- 
ly spaced logarithmically and we interpolate the power 
law from there. The distributions for stream lengths and 
areas must therefore have structures that when combined 
across orders produce smooth power laws. 



For the example of the stream segment length l t 

(s) 



Horton's laws state that the grows by a factor 

of R^s) with each integer step in order lo. In considering 

P(l u , u>), the underlying probability distribution func- 

(s) 

tion for l<s , we postulate that Horton's laws apply for 
every moment of the distribution and not just the mean. 
This generalization of Horton's laws may be encapsulated 
in a statement about the distribution P(iJ~ s \lo) as 

P(lJ s \io) = c l{s) (R n R l{s) )--F l{s) (lJ^R^). (4) 

The factor of (-Rn) - " indicates the that 
J^) =0 d/( S )P(/J s) ,w) cx (Rn)-", i.e., the frequency 
of stream segments of order lo decays according to 
Horton's law of stream number given in equation (Q). 
Similarly, for l u , and we write 



P(L,w) = ci{R n R l{s) ] 



and 



P(a u ,u) = c a (R 2 n )-"F a (a u R- u ), 



P{nn,u) = CniRnf-^F^nn^Rn"), 



(5) 



(6) 



(7) 



where constants c,( s ), q, c a and c„ are appropriate nor- 
malizations. We have used the subscripted versions of the 
lengths and areas, I J- ,l u , and a u , to reinforce that these 
parameters are for points at the outlets of order lo basins 
only. The quantity Uq >uj is the number of streams of order 
lo within a basin of order Ct. This will help with some 
notational issues later on. The form of the distribution 



functions F,(a) > Fi, F a and F n and their interrelationships 
become the focus of our investigations. Since scaling is 
inherent in each of these postulated generalizations of 
Horton's laws, we will often refer to these distribution 
functions as scaling functions. 

We further postulate that distributions of stream seg- 
ment lengths are best approximated by exponential dis- 
tributions. Empirical evidence for this will be provided 
later on in Section M. The normalized scaling function 
F ; ( s ) (u) of equation (^) then has the form 



F lW («) = ie-»/«=F lW (u;0, 



(8) 



where we have introduced a new length scale £ and stated 
its appearance with the notation i*l( s ) («; £). The value of 
£ is potentially network dependent. As we will show, dis- 
tributions of main stream lengths, areas and stream num- 
ber are all dependent on £ and this is the only additional 
parameter necessary for their description. Note that £ is 
both the mean and standard deviation of F ^ (it; £), i.e., 
for exponential distributions, fluctuations of a variable 
are on the order of its mean value. We may therefore 
think of £ as a fluctuation length scale. Note that the 
presence of exponential distributions indicates a random- 
ness in the physical distribution of streams themselves 
and this is largely the topic of our third paper pE . 
Since main stream lengths are combinations of stream 

segment lengths, i.e. l^ = Yli=i lo> i we have that the 
distributions of main stream lengths of order lo basins 
are approximated by convolutions of the stream segment 
length distributions. For this step, it is more appropriate 
to use conditional probabilities such as P(lu\u>) where 
the basin order lo is taken to be fixed. We thus write 

P(L\lo) = P(&K\1) * P(/ (s) 2 |2) * • ■ •PflJ-V)- (9) 
where * denotes convolution. Details of the form 



obtained are given in Appendix A 1 



The next step takes us to the power law distribution 
for main stream lengths. Summing over all stream orders 
and integrating over u — l u we have 



°° />OC 

P(l) ~ ^ / AuP{u,u) 



(10) 



where we have returned to the joint probability for this 
calculation. The integral over u is replaced by a sum 
when networks are considered on discrete lattices. Note 
that the probability of finding a main stream of length 
I is independent of any sort of stream ordering since it 
is defined on an unordered network. The details of this 



calculation may be found in Appendix A 2 where it is 
shown that a power law P(l) cx follows from the 
deduced form of the P{l u ,uS) with 7 = In R n / lni? ( s ) . 

IV. ESTIMATION OF HORTON RATIOS 

We now examine the usual Horton's laws in order to 
estimate the Horton ratios. These ratios are seen as 
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FIG. 2: Horton's laws for the order £1 = 11 Mississippi river 
basin network. For (a), the ordinate axis is logarithmic (base 
10) representing number for stream number n„ (circles), km 2 
for area a u (squares), and km for both main stream length 

- fs) 

Iuj (triangles) and stream segment length l u (diamonds). 

Note the good agreement between l u and I J- . In (b), the 
stream number data (circles) has been inverted from that 
in (a), i.e., the plot is of nj 1 . This is compared with the 
dimensionless a^/aa (squares) showing good support for the 
prediction the slopes are equal, i.e., R a = Rn- 



intrinsic parameters in the probability distribution func- 
tions given above in equations (f|), (g), (||) and (Q). 

Figure |](a) shows the stream order averages of ZW, Z, 
a and n for the Mississippi basin. Deviations from expo- 
nential trends of Horton's laws are evident and indicated 
by deviations from straight lines on the semi-logarithmic 
axis. Such deviations are to be expected for the smallest 
and largest orders within a basin ||, |l). For the small- 
est orders, the scale of the grid used becomes an issue 
but even with infinite resolution, the scaling of lengths, 
areas and number for low orders cannot all hold at the 
same time For large orders, the decrease in sample 
space contributes to these fluctuations since the number 
of samples of order u> streams decays exponentially with 
order as {Rnf 1 ^ 1 ^ . Furthermore, correlations with over- 
all basin shape provide another source of deviations [pi] . 
Nevertheless, in our theoretical investigations below we 
will presume exact scaling. Note also that the equiva- 
lence of R n and R a is supported by Figure ||(b) where 
the stream numbers n w have been inverted for compari- 
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1 ABLE I: Hoggin ratios for the Mississippi River []36| 

tch range ofRorders (u!i,u>2), estimates of the ratios 
o jtained mm/ simple regression analysis. For each 
t; , a meafti, standard deviation a and normalized devl 
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these estimates but not all are shown The values 
ft r Ri are especially robust while some variation is obse 
'or t 



betwi 



estunatcs ot 
n trre raMc 



„ and ha- Uoo 
ana i?Pand'als< 



For 
are 
qu^nti- 
tion 
ised 
obt a|ined 
ved 



pf merit is observed 
weerrff; ana R r s \ . 









CO 








u) range 


Rn 


Ra 


Ri 


%) 


Ra/ Rn 


Ri/R l(l 


[2,3] 


5.05 


4.69 


2.10 


1.65 


0.93 


1.28 


[2,5] 


4.65 


4.64 


2.11 


1.92 


1.00 


1.10 


[2,7] 


4.54 


4.63 


2.16 


2.11 


1.02 


1.03 


[3,4] 


4.54 


4.73 


2.10 


2.01 


1.04 


1.05 


[3,6] 


4.51 


4.62 


2.15 


2.15 


1.02 


1.00 


[3,8] 


4.44 


4.55 


2.19 


2.23 


1.02 


0.98 


[4,6] 


4.52 


4.59 


2.18 


2.24 


1.02 


0.97 


[4,8] 


4.42 


4.51 


2.21 


2.27 


1.02 


0.97 


[5,7] 


4.39 


4.62 


2.25 


2.39 


1.05 


0.94 


[6,7] 


4.19 


4.55 


2.26 


2.40 


1.09 


0.94 


[7,8] 


4.50 


4.21 


2.15 


2.12 


0.94 


1.02 


mean fi 


4.51 


4.58 


2.17 


2.15 


1.01 


1.02 


std dev a 


0.17 


0.12 


0.05 


0.19 


0.03 


0.08 




0.038 


0.026 


0.024 


0.089 


0.034 


0.078 



TABLE II: 
per Table f 



Horton ratios for the Amazon 



Details are as 



son. Similar agreement is found for the Amazon and Nile 
as shown in Tables |, |l[ and II] which we now discuss. 

Table [j] shows the results of regression on the Missis- 
sippi data for various ranges of stream orders for stream 
number, area and lengths. Tables |l| and [II show the 
same results carried out for the Amazon and Nile. Each 
table presents estimates of the four ratios R n , R a , Ri 
and R^s) ■ Also included are the comparisons R a /R n and 
i?z/i? z ( s ), both of which we expect to be close to unity. 
For each quantity, we calculate the mean //, standard 
deviation a and normalized deviation a / '/i. 

Note the variation of exponents with choice of order 
range. This is the largest source of error in the calcula- 
tion of the Horton ratios. Therefore, rather than taking a 
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TABLE III: Horton ratios for the Nile 
per Table | Here 2 < u>i < oj 2 < 7. 



Details are as 



single range of stream orders for the regression, we exam- 
ine a collection of ranges. Also, the deviations for high 
and low orders observed in Figures |(a) and |(b) do of 
course affect measurements of the Horton ratios. In all 
cases, we have avoided using data for the smallest and 
largest orders. 

For the three example networks given here, the state- 
ments R a = R n and Ri = R^s) are well supported. The 
majority of ranges give R n /R a and Ri/R^ very close to 
unity. The averages are also close to one and are different 
from unity mostly by within 1.0 and uniformly by within 
1.5 standard deviations. 

The normalized deviations, ie., cr/yU, for the four ratios 
are all below 0.05. No systematic ordering of the a / '/i is 
observed. Of all the data, the values for Ri in the case 
of the Mississippi are the most notably uniform having 
a I [i = 0.015. Throughout there is a slight trend for 
regression on lower orders to overestimate and on higher 
orders to underestimate the average ratios, while reason- 
able consistency is found at intermediate orders. 

Thus, overall the ranges chosen in the tables give a 
reasonably even set of estimates of the Horton ratios and 
we will use these averages as our estimates of the ratios. 



V. EMPIRICAL EVIDENCE FOR HORTON 
DISTRIBUTIONS 

A. Stream segment length distributions 

We now present Horton distributions for the Mississip- 
pi, Amazon, and Nile river basins as well as the Schei- 
degger model. Scheidegger networks may be thought 
of as collections of random-walker streams and are ful- 
ly defined in |l(| and extensively studied in ^l). The 
forms of all distributions are observed to be the same in 
the real data and in the model. 

The first distribution is shown in Figure ||(a). This 
is the probability density function of 1^4,, fourth order 
stream segment lengths, for the Mississippi River. Dis- 
tributions for different orders can be rescaled to show 
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FIG. 3: Plot (a) shows an example distribution of stream seg- 
ment lengths, P(lJ- b \uj), for the Mississippi for order uj = 4. 
The lengths here are in kilometers. The semi-logarithmic axes 
indicate the distribution is well approximated by an expo- 
nential. The value of the length scale £ (see equation (^)) is 
estimated to be approximately 800 meters. Rescaled versions 
of the same stream segment length distributions for uj — 3 
(circles), uj — 4 (squares), uj = 5 (triangles), and uj — 6 (dia- 
monds), are shown in (b). The rescaling is done according 
to equation (^). The values of the Horton ratios used are 
R n = 4.69 and R,( s ) = 2.33 as determined from Table f| 



satisfactory agreement. This is done using the postulat- 
ed Horton distribution of stream segment lengths given 
in equation ([|). The rescaling is shown in Figure ||(b) 
and is for orders uj = 3, . . . , 6. Note the effect of the 
exponential decrease in number of samples with order 
is evident for uj = 6 since P(l^ s \) is considerably scat- 
tered. Nevertheless, the figure shows the form of these 
distributions to be most closely approximated by expo- 
nentials. We observe similar exponential distributions for 
the Amazon, the Nile and the Scheidegger model. The 
fluctuation length scale £ is found to be approximately 
800 meters for the Mississippi, 1600 meters for the Ama- 
zon and 1200 meters for the Nile. 

Since £ is based on the definition of stream ordering, 
comparisons of £ are only sensible for networks that are 
measured on topographies with the same resolution. The 
above values of £ are approximate and our confidence 
in them would be improved with higher resolution data. 
Nevertheless they do suggest that fluctuations in net- 
work structure increase as we move from the Mississippi 
through to the Nile and then the Amazon. 
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FIG. 4: Plot (a) shows an example distribution for order 
ui = 5 main stream lengths (measured in km) for the Ama- 
zon. The distribution is unimodal with what is a reasonable 
approximation of an exponential tail. In (b), distributions of 
main stream length for for uj = 3 (circles), u) = 4 (squares), 
uj — 5 (triangles), and uj = 6 (diamonds), are rescaled accord- 
ing to equation (ph . The values of the Horton ratios used here 
are R n = 4.51 and -R ; ( s ) — 2.17, taken from Table 



FIG. 5: Given in (a) is an example distribution of order 
ui — 6 main stream lengths for the Scheidegger model. The 
same form is observed as for real networks such as the Amazon 
(Figure^). In the same way as Figure [|(b), (b) show rescaled 
distributions of main stream length for for u> = 4 (circles) , u) — 
5 (squares), ui = 6 (triangles), and uj — 7 (diamonds). Note 
that in (b), distributions are not normalized with respect to a 
fixed basin order £1 and hence the vertical offset is arbitrary. 
The values of the ratios used here are R„ ~ 5.20 and Ri ~ 
3.00 H. 



B. Main stream segment length distributions 

The distributions of uj = 4 main stream lengths for 
the Amazon River is shown in Figure |^(a). Since main 
stream lengths are sums of stream segment lengths, their 
distribution has a single peak away from the origin. How- 
ever, these distributions will not tend towards a Gaus- 
sian because the individual stream length distributions 
do not satisfy the requirements of the central limit theo- 
rem |3^] . This is because the moments of the stream seg- 
ment length distributions grow exponentially with stream 
order. As the semi-logarithmic axes indicate, the tail may 
be reasonably well (but not exactly) modeled by expo- 
nentials. There is some variation in the distribution tails 
from region to region. For example, corresponding dis- 
tributions for the Mississippi data do exhibit tails that 
are closer to exponentials. However, for the present work 
where we are attempting to characterize the basic forms 
of the Horton distributions, we consider these deviations 
to be of a higher order nature and belonging to the realm 
of further research. 

In accord with equation (||), Figure [|(b) shows the 
rescaling of the main stream length distributions for 
uj = 3, . . . , 6. The ratios used, R n — 4.49 and Ri = 
2.19(~ R (s) = 2.17) are taken from Table || Given 



the scatter of the distributions, it is unreasonable to per- 
form minimization techniques on the rescaled data itself 
in order to estimate R n and Ri. This is best done by 
examining means, as we have done, and higher order 
moments which we discuss below. Furthermore, varying 
R n and Ri from the above values by, say, ±0.05 does not 
greatly distort the visual quality of the "data collapse." 

Similar results for the Scheidegger model are shown in 
Figure ||. The Scheidegger model may be thought of as 
a network defined on a triangular lattice where at each 
lattice site one of two directions is chosen as the stream 
path |ll| ^l| . Figure |^(a) gives a single example distribu- 
tion for main stream lengths of order uj = 6 basins. The 
tail is exponential as per the real world data. Figure ||(b) 
shows a collapse of main stream length distributions for 
orders uj = 4 through 7. In contrast to the real data 
where an overall basin order is fixed (f2), there is no max- 
imum basin order here. The distributions in Figure ||(b) 
have an arbitrary normalization meaning the absolute 
values of the ordinate are also arbitrary. Otherwise, this 
is the same collapse as given in equation (||). For the 
Scheidegger model, our simulations yield R n ~ 5.20 and 
i? ; ( s ) ~ 3.00 For all distributions, we observe similar 
functional forms for real networks and the Scheidegger 
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FIG. 6: The distribution of drainage areas for lu — 4 sub- 
basins of the Nile are shown in (a) . All areas are measured in 
km 2 . An exponential tail is observed as per the distributions 
of stream segment length (Figure 0) and main stream length 
(Figure ^). In (b), distributions or drainage area for for lo — 
3 (circles), lo = 4 (squares), lo — 5 (triangles), and lo — 
6 (diamonds), are rescaled according to equation (^K The 
rescaling uses the estimate R n = 4.42 found in Table III. 



model, the only difference lying in parameters such as 
the Horton ratios. 



FIG. 7: Summation of main stream length distributions for 
the Mississippi. Both axes are logarithmic,the unit of length 
is km and the vertical axis is probability density with units 
of km -1 . Distributions of for orders lo = 3 (circles), lo = 4 
(squares), and to = 5 (triangles), are shown. As expected, the 
distributions sum together to give a power law tail (stars). 
The power law distribution (which is vertically offset by an 
order of magnitude for clarity) is the summation of the dis- 
tributions below as well as the distribution for order lo — 6 
main stream lengths. 



C. Drainage area distributions 



D. Summing distributions to form power laws 



Figure |g shows more Horton distributions, this time for 
drainage area as calculated for the Nile river basin. In 
Figure an example distribution for to = 4 sub-basins is 
presented. The distribution is similar in form to those of 
main stream lengths of Figure f|, again showing a reason- 
ably clear exponential tail. Rescaled drainage area dis- 
tributions for lv = 3, . . . , 6 are presented in Figure ||(b). 
The rescaling now follows equation @ . Note that if R n 
and R a were not equivalent, the rescaling would be of the 
form 



(11) 



Since we have asserted that R n = R a , equation fll] 
reduces to equation @. The Horton ratio used here is 
R n = 4.42 which is in good agreement with R a = 4.53, 
the respective standard deviations being 0.17 and 0.10. 



As stated in Section [II, the Horton distributions of 
and l u must combine to form power law distributions 
for a and I (see equations || and ]Io| ) . Figure [?] provides 
empirical support for this observation for the example 
main stream lengths of the Mississippi network. The 
distributions for lo = 3, 4 and 5 main stream lengths 
are individually shown. Their combination together with 
the distribution of Iq gives the reasonable approximation 
of a power law as shown. The area distributions com- 
bine in the same way. Note that the distributions do not 
greatly overlap. Each point of the power law is therefore 
the addition of significant contributions from only two or 
three of the separate distributions. The challenge here 
then is to understand how rescaled versions of Fi, being 
the basic form of the P(l u ,u)), fit together in such a clean 
fashion. The details of this connection are established in 



Both figures are taken from the data of Table III 



Appendix A 2 
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E. Connecting distributions of number and area 

In considering the generalized Horton distributions for 
number and area, we observe two main points: a calcula- 
tion in the vein of what we are able to do for main stream 
lengths is difficult; and, the Horton distributions for area 
and number are equivalent. 

In principle, Horton area distributions may be derived 
from stream segment length distributions. This follows 
from an assumption of statistically uniform drainage den- 
sity which means that the typical drainage area drained 
per unit length of any stream is invariant in space. Apart 
from the possibility of changing with space which we will 
preclude by assumption, drainage density does naturally 
fluctuate as well O], Thus, we can write a ~ pJ2u 
where the sum is over all orders and all stream segments 
and p is the average drainage density. 

However, wc need to know for an example basin, how 
many instances of each stream segment occur as a func- 
tion of order. For example, the number of first order 
streams in an order fl basins is nn,x- Given the distri- 
bution of this number, we can then calculate the distri- 
bution of the total contribution of drainage area due to 
first order streams. But the distributions of hq^ are not 
independent so we cannot proceed in this direction. 

We could potentially use the typical number of order 
ui streams, (i? n ) n_w . Then the distribution of total area 
drained due to order u> streams would approach Gaussian 
because the individual distribution are identical and the 
central limit theorem would apply. However, because the 
fluctuations in total number of stream segments are so 
great, we lose too much information with this approach. 
Indeed, the distribution of area drained by order uj stream 
segments in a basin reflects variations in their number 
rather than length. Again, we meet up with the prob- 
lem of the numbers of distinct orders of stream segment 
lengths being dependent. 

One final way would be to use Tokunaga's law g, f«], 
EjJ [Is). Tokunaga's law states that the number of 
order v side branches along an (absorbing) stream seg- 
ment of order /1 is given by 

T fc = T 1 (i2 lW ) fc + 1 . (12) 

where k = fj, — v. The parameter T% is the average num- 
ber of side streams having order v = p, — 1 for every order 
fi absorbing stream. This gives a picture of how a net- 
work fits together and may be seen to be equivalent to 
Horton's laws ||. Now, even though we also understand 
the distributions underlying Tokunaga's law ||l6| , similar 
technical problems arise. On descending into a network, 
we find the number of stream segments at each level to 
be dependent on all of the above. 

Nevertheless, we can understand the relationship 
between the distributions for area and number. What 
follows is a generalization of the finding that R n = R a . 
The postulated forms for these distributions were giv- 
en in equations @ and (]?]). Consider nn,i, the num- 




12 3 4 

x 10 4 

FIG. 8: Comparison of number and area distributions for the 
Scheidegger model. Area is in terms of lattice units. In the 
inset plot, the raw distributions shown are P(aa\6) (circles) 
and P(ri6,i|6) (continuous line). The latter is the probability 
of finding ne,i source streams in an order lu = 6 basin. In the 
main plot, the number distribution has been rescaled to be 
l/4P(ne,i|6) as a function of 4ri6,i and the area distribution 
is unrescaled (the symbols are the same as for the inset plot). 
For the Scheidegger model, source streams occur at any site 
with probability of 1/4, hence the rescaling by a factor of four. 



ber of first order streams in an order f2 basin. Assum- 
ing that, on average, first order streams are distributed 
evenly throughout a network, then this number is simply 
proportional to an- As an example, Figure ^ shows data 
obtained for the Scheidegger model. For the Scheidegger 
model, first order streams are initiated with a 1/4 proba- 
bility when the flow at the two upstream sites is randomly 
directed away, each with probability 1/2. Thus, for an 
area an, we expect and find riQ tU , = an/4. 

For higher internal orders, we can apply a simple 
renormalization. Assuming a system with exact scaling, 
the number of streams hq^ is statistically equivalent to 
nn-u)+i,i- Since the latter is proportional to an-w+i we 
have that 

(13) 

where the constant of proportionality is the density of 
order u> streams, Clearly, this equivalence improves as 
number increases, i.e., the difference O — ui increases. 

While we do not have exact forms for the area or num- 
ber distributions, we note that they are similar to the 
main stream length distributions. Since source streams 
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Mississippi: Horton moments 




1 2 3 4 5 6 7 8 
order co 



10 11 



FIG. 9: A comparison of moments calculated for main stream 
length distributions for the Mississippi River. 



are linear basins with the width of a grid cell, the dis- 
tribution of ai is the same as the distribution of l\ and 
l^ s \, a pure exponential. Hence, nn t n-i is also an expo- 
nential. For increasing uj, the distribution of a w becomes 
single peaked with an exponential tail, qualitatively the 
same as the main stream length distributions. 



!(/ (s) i) 



Since the characteristic length of order ui 



streams is (i? ; ( s ))" , we therefore have 

((L {s) ) q ) = M q (R l ( S )) { "- 1)q = ?! (L (s) ) q ■ (15) 

Since main stream lengths are sums of stream segment 
lengths, so are their respective moments. Hence, 



- £((/ (s) *)«), 

fc=l 



fc=l 



k=l 



R 



1 



(16) 



We can now determine the log-space separation of 
moments of main stream length. Using Stirling's approx- 
imation p4| that Inn! ~ {n + 1/2) Inn — n we have 



ln((U 9 )~5 € + (J2. W r + In g 



(17) 



where C is a constant. The lng term inside the square 
brackets in equation [l?] creates small deviations from lin- 
earity for 1 < uj < 15. Thus, in agreement with Figure]^, 
we expect approximately linear growth of moments in 
log-space. 



VI. HIGHER ORDER MOMENTS 



VII. LIMITATIONS ON THE PREDICTIVE 
POWER OF HORTON'S LAWS 



Finally, we discuss the higher order moments for the 
generalized Horton distributions. Figure || presents 
moments for distributions of main stream lengths for the 
case of the Mississippi. These moments are calculated 
directly from the main stream length distributions. A 
regular logarithmic spacing is apparent in moments for 
orders ranging from 3 to 7. 

To see whether or not this is expected, we detail a 
few small calculations concerning moments starting from 
the exponential form of stream segment lengths given in 
equation (||). As noted previously, for an exponential 
distribution, F ; ( s ) (u) = l- l e- u /t, the mean is simply 
(u) = £. In general, the gth moment of an exponential 
distribution is 



00 mi 



u=0 £ 



£, q / x q e~ x dx = ql4 q . (14) 

Jx=0 



Assuming scaling holds exactly for across all orders, 
the above is precisely ((/ (s) i) 9 ). Note that ((^ (s) i) 9 ) = 



In this last section, we briefly examine deviations from 
scaling within this generalized picture of Horton's laws. 
The basic question is given an approximate scaling for 
quantities measured at intermediate stream orders, what 
can we say about the features of the overall basin? 

As noted in the previous section, all moments of the 
generalized Horton distributions grow exponentially with 
order. Coupling this with the fact that ex R~ u , i.e., 
the number of samples of order uj basins decreases expo- 
nentially with uj, we observe that a basin's a and I will 
potentially differ greatly from values predicted by Hor- 
ton's laws. 

To illustrate this, Figure [l(] specifically shows the dis- 
tributions P(h) and P(h) scaled up to give P(hi) for 
the Congo river. The actual Congo's length measured 
at this 1000 meter resolution is represented by the sol- 
id line and is around 57% of the distribution's mean as 
indicated by the dashed line. Nevertheless, we see that 
the measured length is within a standard deviation of the 
predicted value. 

we provide a comparison of predicted 



In Table [V 



versus measured main stream lengths and areas for the 



10 



x 10 




/ R, (m) 

CO / 
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FIG. 10: Comparing the generalized Horton length distri- 
bution rescaled to the level of order Q — 11 basins with the 
Congo River itself. The two distributions are for orders u — 3 
(squares) and u = 4 (circles) stream lengths and the Horton 
ratio is estimated to be Ri = 2.39 j45|. The dashed line rep- 
resents the mean of these scaled up distributions while the 
solid line marks In, the measured length of the Congo at a 
1000 meter resolution. The actual length is within a standard 
deviation of the mean being around 50% of In. Table |Fv| 
shows comparisons for various river networks for both area 
and length data. 



basin: 



Mississippi 

Amazon 

Nile 

Congo 

Kansas 



In In <Tj In /In oi/ln 



4.92 11.10 5.60 0.44 0.51 

5.75 9.18 6.85 0.63 0.75 

6.49 2.66 2.20 2.44 0.83 

5.07 10.13 5.75 0.50 0.57 

1.07 2.37 1.74 0.45 0.73 



an an 



an /an 



2.74 7.55 5.58 0.36 

5.40 9.07 8.04 0.60 

3.08 0.96 0.79 3.19 

3.70 10.09 8.28 0.37 

0.14 0.49 0.42 0.28 



of stream lengths and areas. 

Thus, we see that by using a probabilistic point of view, 
this generalized notion of Horton's laws provides a way of 
discerning the strength of deviations about the expected 
mean. In general, stronger deviations would imply that 
geologic conditions play a more significant role in setting 
the structure of the network. 



VIII. CONCLUSION 

The objective of this work has been to explore 
the underlying distributions of river network quantities 
defined with stream ordering. We have shown that func- 
tional relationships generalize all cases of Horton's laws. 
We have identified the basic forms of the distributions for 
stream segment lengths (exponential) and main stream 
lengths (convolutions of exponentials) and shown a link 
between number and area distributions. Data from the 
continent-scale networks of the Mississippi, Amazon, and 
Nile river basins as well as from Scheidegger's model of 
directed random networks provide both agreement with 
and inspiration for the generalizations of Horton's laws. 
Finally, we have identified a fluctuation length scale £ 
which is a reinterpretation of what was previously iden- 
tified as only a mean value. We see the study of the 
generalized Horton distributions as integral to increas- 
ing our understanding of river network structure. We 
also suggest that practical network analysis be extended 
to measurements of distributions and the length scale £ 
with the aim of refining our ability to distinguish and 
compare network structure. 

By taking account of fluctuations inherent in network 
scaling laws, we are able to see how measuring Horton's 
laws on low-order networks is unavoidably problematic. 
'Moreo ver, as we have observed, the measurement of the 



rtiJAon ratios is in general a delicate operation suggesting 
t9ii$ many previous measurements are not without error. 
0-¥he theoretical understanding of the growth and evolu- 
tftS? of river networks requires a more thorough approach 
t@i8measurement and a concurrent improvement in the 



TABLE IV: Comparison of predicted versus measured main 
stream lengths for large scale river networks. The dimensions 
of all lengths and areas are 10 6 m and 10 12 m 2 respectively. 
Here, In, is the actual main stream length of the basin, In 
the predicted mean value of In , <Ji the predicted variance and 
ci/ln the normalized deviation. The entries for the basin area 
data have corresponding definitions. 



statistical description of river network geometry. The 
present consideration of a generalization of Horton's laws 
is a necessary step in this process giving rise to stronger 
tests of both real and synthetic data. In the following 
paper we round out this expanded picture of net- 
work structure by consdering the spatial distribution of 
network components. 



basins studied here. The mean for the scaled up distribu- 
tions overestimates the actual values in all cases except 
for the Nile. Also, apart from the Nile, all values are 
within a standard deviation of the predicted mean. The 
coefficients of variation, cr a /an and <Ji/ln, all indicate 
that fluctuations are on the order of the expected values 
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APPENDIX A: ANALYTIC CONNECTIONS 
BETWEEN STREAM LENGTH DISTRIBUTIONS 

In this appendix we consider a series of analytic cal- 
culations. These concern the connections between the 

(s) 

distributions of stream segment lengths lj~ , ordered 
basin main stream lengths l u and main stream lengths 
I. We will idealize the problem in places, assuming per- 
fect scaling and infinite networks while making an occa- 
sional salubrious approximation. Also, we will treat the 
problem of lengths fully noting that derivations of dis- 
tributions for areas follow similar but more complicated 
lines. 

We begin by rescaling the form of stream segment 
length distributions 



(R n - l)(R n R l{s) )- u F l{s) (lR^). (Al) 



The normalization e,( s ) 
ment that 



R n — 1 stems from the require- 



^ W («) = l, 



(A2) 



which is made purely for aesthetic purposes. As we have 
suggested in equation (||) and demonstrated empirical 
support for, (u) is well approximated by the expo- 
nential distribution S i ~ 1 e~ u ^. For low u and also we have 
noted that deviations do of course occur but they are suf- 
ficiently insubstantial as to be negligible for a first order 
treatment of the problem. 



Distributions of main stream lengths as a 
function of stream order 



We now derive a form for the distribution of main 
stream lengths P(Z w |o;). As we have discussed, since 



luj = Ya=i lw B \ we have the convolution (^). The right- 
hand side of equation (Q) consists of exponentials as per 
equation (||) so we now consider the function K u (u;a) 
given by 

K u (u; a) = a ie - aiu * a 2 er a * u a^er a " u , (A3) 

where a = (ax, a%, . . . , a u ). We are specifically interested 
in the case when no two of the ai are equal, i.e., a* 7^ aj 
for all i 7^ j. To compute this w-fold convolution, we 
simply examine the K u (u; a) for u> — 2 and u> = 3 and 
identify the emerging pattern. For a — (ai, a 2 ) we have, 
omitting the prefactors for the time being, 



(A4) 



at — a 2 ai — a 2 a 2 — ai 

providing a\ 7^ a 2 . Convolving this with e~ a3U we obtain 

' p — aiu _ p — 02« ' 
e -a lU ^ e -a 2 u ^ e -a 3 u = [ ] ^ e - as n 



a\ — a 2 

p —a 2 u _ p—d3U 



(ai - a 2 )(ai - a 3 ) (ai - a 2 )(a 2 - a 3 )' 



(ai - a 2 )(ai - a 3 ) 



(02 - ai)(a 2 - a 3 ) (a 3 - ai)(a 3 - a 2 ) 
Generalizing from this point, we obtain 



K u {u;a) = rQa, 



\i— 1 / i 



(A5) 



(A6) 



Now, setting ai = l/(£(i? ; ( s )) 1 r ) and carrying out some 
manipulations we obtain the following expression for 



1 1 y 



y- 1 



1 r-^-w^r- 1 ^^^ 



E 



1 if.^wr^wJ^^wrVa.w 



E 



_J L_ V e-'-^^iW r 1 W ( "~ 2) (A7) 



Note that we have added in a factor of l/{R n ) w for the appropriate normalization. In addition, one observes 
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that P(0,ui) = for all uj > 1 since all convolutions of 
pairs of exponentials vanish at the origin. Furthermore, 
the tail of the distribution is dominated by the exponen- 
tial corresponding to the largest stream segment. 

The next step is to connect to the power law distribu- 
tion of main stream lengths, P(l) (see Figure ^ and the 
accompanying discussion). On considering equation (|Io| ) 
we see that the problem can possibly be addressed with 
some form of asymptotic analysis. 

Before attacking this calculation however, we will sim- 
plify the notation keeping only the important details of 
the P(L), u>). Our main interest is to see how equa- 
tion (nq) gives rise to a power law. We transform the 
outcome of equation (A7) by using n — u>, u — 



i? ; ( s ), and s = R n , neglecting multiplicative con- 
stants and then summing over stream orders to obtain 



1 ^ r (n-2)i e — u/r' 1 



(A8) 



The integration over l u has been omitted meaning that 
the result will be a power law with one power lower than 
expected. 



2. Power law distributions of main stream lengths 

We no w sh ow that this sum of exponentials G{u) in 
equation ( |A8| ) does in fact asymptotically tend to a power 
law. We first interchange the order of summation replac- 
m S E„=i E i= i Wlth E i= i E„=* to give 

r {n— 2)i 



G{u) = Ve'^'T * 



(A9) 



i=i 



We thus simply have a sum of exponentials to contend 
with. The coefficients Cj appear unwieldy at first but 
do yield a simple expression after some algebra which we 
now perform: 



,.(71—2)1 



^s"n;=i,^-r*) 

1 ~ r {n-2)i 



r li 2), » g i f A„ 2)i r -I/' -2), 



Z-^l qTl- 



r (n—i)i 



1 1 1 




A10) 



In reaching the last line we have shifted the indices in 
several places. In the last bracketed term we have set 
k = j — i and then m = n — i while in the first bracketed 
term, we have used —k = j — i. Immediately of note is 
that the last term is independent of i and may thus be 
ignored. 

The first bracketed term does depend on i but con- 

') we have 



verges rapidly. Writing Di = nl=i(l — r h '■ 

that Di — D m Ilfc=m (1 — r ~ fe ) • Taking m to be fixed and 

"~ k is approximated well by 
exp{— r~ k } for k > to, we then have 



large enough such that 1 



Di 



D m exp 



D n 



exp 




(All) 



As i — > oo, Di clearly approaches a product of D m and 
a constant. Therefore, the first bracketed term in equa- 
tion (A10) may also be neglected in an asymptotic anal- 
ysis. 

Hence, as i — > oo, the coefficients Ci are simply given 

by 



Ci oc 



1 



(A12) 



and we can approximate G(u) as, boldly using the equal- 
ity sign, 



G(u) = AS(u) 



A 



00 g-u/r' 

X] r i(l+7) ' 
i=0 



(A13) 



where A comprises the constant part of the Ci and factors 
picked up by shifting the lower limit of the index i from 
1 to 0. We have also used here the identification 

s = r 7 . (A14) 

We turn now to the asymptotic behavior of S(u), this 
being the final stretch of our analysis 

There are several directions one may take at this 
point. We will proceed by employing a transformation 
of S(u) that is sometimes referred to as the Sommerfeld- 
Watson transformation and also as Watson's lemma J4t| , 
p. 239]. Given a sum over any set of integers /, say 
S = Ene//( n )' it can be written as the following inte- 
gral 



S = 



1 

2tH 



IT COS 7TZ 

sin7rz 



f(z)dz. 



(A15) 
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where C is a contour that contains the points on the real 
axis n + iO where n £ I and none of the points of the 
same form with n £ Z/7. Calculation of the residues of 
the simple poles of the integrand return us to the original 
sum. 

Applying the transformation to S(u) we obtain 

S (u) = — I I^Ile-^-'r-'^dz. (A16) 
2iri J c simrz 

The contour C is represented in Figure 




FIG. 11: Co ntou r C used for evaluation of the integral given 
in equation (A16). The poles are situated at n + Oi where 
"£{0,1,2,4 



We hrst make a change of variables, r~ z — p. Substi- 
tuting this and dz = — dp/phir into equation ( A16| ) we 
have 



S(u) 



1 

2ttI 
1 



' 7r cos — 7r In pi lnr 



sin — 7r In pj lnr 
' 7TCOS7rlnp/lnr 



e - u Pp( 1+ i\-dp/plnr) 



-up 7 



p~<dp. 



(A17) 



2ilnr J c sin 7r Inp/ lnr 

The transformed contour C is depicted in Figure [l^. 

As u — > oo, the contribution to integral from the neigh- 
borhood of p = dominates. The introduction of the sin 
and cos terms has created an interesting oscillation that 
has to be handled with with some care. We now deform 
the integration contour C into the contour C" of Fig- 
ure [13 focusing on the interval along the imaginary axis 
[—*,!]. Choosing this path will simplify the cos and sin 
expressions which at present have logs in their arguments. 

The integral S(u) is now given by S(u) ~ I(u) + c.c. 
where 



I(u) 



-1 



7TCOS7rlnp/lnr _ u 



e- up p' f dp. (A18) 



2ilnr J sin7rln/5/lnr 

Writing p = a + ir with a = 0, we have dp = idr and 
the following for the cos and sin terms: 



cos 7r In pj lnr 



pi-it j In r _j_ p^i^ / In r 



r iir J \nr p — 7r /21nr 



+ T 



-27T / In r p 7T /2 In r 



(A19) 




FIG. 12: C ontou r C' used for evaluation of the integral given 
in equation (A17) as deduced from contour C (Figure p"l|) with 
the transformation p = r~ z . The negative real axis is a branch 
cut. 




FIG. 13: Contou r C" used for evaluation of the integral 
given in equation (A16). The poles are situated at n + Oi 
where n £ {0, 1,2,.. .}. 



and 



sin7rlnp/ lnr = 



pin / In r p— vk j In ? 

2i 



j-iir/ In r 7r 2 /2 In r — in/ In r ^tz' 2 /2 In r 

2i ' 

The cot term in the integrand becomes 

COSTThWlnr 1 + T ^/^r e --K 2 /\nr 

-i- 



(A20) 



sin7rln pj lnr 



' ^ j-2in j In rg— 7r 2 / In r 

1 + S(t) 

•i-s(ry 



(A21) 



where *(r) = T 2^/inr e -7r 2 /inr The integral I(u) now 
becomes 



I(u) 



i / >1 1 + <5(t) 



21nr J 1 - S(t) 



-iiiT T 7 e i7T7/2 ( j T 
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I e -^ r ^ ! + 1 dr. (A22) 



2 lnr 



IUT_7 " 

l-<5(r)' 



Now, since |<5(r)| = e -* 2 / lnr < 10" 4 (taking r = i? ;(s) 
2.5), we can expand the expression as follows 



1 + 5 
1 - S 



= (l + S)(l + S + S 2 + ...) 
= 1 + 25 + 2S 2 + 25 a + ... 



(A23) 



The integral in turn becomes 



/(«) = 



,■1+7 /■! 



21nr 



^ _|_ 2 T 2i7r/ lni-g-TT 2 / liir _|_ 2 T 4i7r / lnI "g-27r 2 /lnr _|_ 
+ 2r 2m 7 r/lnr e -n7 r 2 /lnr + ^ ( A24 ) 

The basic n-th integral in this expansion is 

I n { U )= f T 7+2«Wlnr e - JU r d ^ 

Jo 

Substituting ut = w and replacing the upper limit w — u 
with w = oo we have 

/>0O 

/„(u) = u -(l+7+2ni7T/lnr) / dW +2m,r/ ln '"e"™' , 



(ra) 
(ra) 



- ( l+7+2?it7r / In r) 



-(l+7+2ni7r/ lnr) 



7+2m7r/ lnr — iw 



idw(iw) 



dv(v)~< +2m7v/lar e- v , 



(m)^( 1+7+2ra7r / lnr )r(7 + 2nin/\nr). 



(A26) 



Here, we have rotated the contour along the imaginary 
iiu-axis to the real v-axis and identified the integral with 
the gamma function T Ej . The integral can now be 
expressed as 



I(u) 



21nru 1+ T 



1 + 2^ u ~ 2nl7T / ln T(7 + 2niir/ lnr) 

n=l 

(A27) 



We now need to show that the higher order terms are 
negligible. Note that their magnitudes do no vanish 
with increasing u but instead are highly oscillatory terms. 
Using the asymptotic form of the Gamma function J17| 



r(» = z 2 - 1 / 2 e- 2 V2^(l + 0(l/z)) . 



we can estimate as follows for large n that 



(A28) 



|r(l + 7 + 2rwr/lnr)| 
- |(2 l7 rn/lnr + 1 + 7 ) 2 ™"/ 1.^+1/2+7^7-1^ 

= I (e l7r / 2 27rn/ lnr) 2 "™/ lnr+1 / 2+ ~< e^^V^l 

= e- 7r2 "/ lnr n 7+1 / 2 (27r/e) 1+ ''(lnr)- 1 / 2 - 7 . (A29) 
Hence, T(l + 7 + 2nin / lnr) vanishes exponentially with 
n. For the first few values of n taking 7 = 3/2 and 
r = 2.5, we have T(l + 7 + 2i7r/lnr) ~ 2.5 x 10~ 3 and 
T(l + 7 + 4i7r/lnr) ~ 2.1 x 1CT 6 showing that these 
corrections are negligible. 

Hence we are able estimate S(u) to first order as 



S(u) ~ — , 
lnr 



,-1-7 



(A30) 



Thus we have determined that a power law follows from 
the initial assumption that stream segment lengths follow 
exponential distributions. 

This equivalence has been drawn as an asymptotic one, 
albeit one where convergences have been shown to be 
rapid. The calculation is clearly not the entire picture 
as the solution does contain small rapidly-oscillating cor- 
rections that do not vanish with increasing argument. A 
possible remaining problem and one for further investi- 
gation is to understand how the distributions for main 
stream lengths l u fit together over a range that is not to 
be considered asymptotic. Nevertheless, the preceding 
is one attempt at demonstrating this rather intriguing 
breakup of a smooth power law into a discrete family of 
functions built up from one fundamental scaling function. 
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